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We consider deeply virtual Compton scattering and deep inelastic scattering in presence of Regge 
exchanges that are part of the non-perturbative quark-nucleon amplitude. In particular we discuss 
contribution from the Pomeron exchange and demonstrate how it leads to Regge scaling of the 
Compton amplitude. Comparison with HERA data is given. 
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I. INTRODUCTION 

In the past two decades, notable theoretical activity 
has been dedicated to the study of the generalized parton 
distributions (GPD's) [THE]- GPD's allow one to access 
the nucleon structure in a more detailed manner than 
the parton distribution functions (PDF's) studied within 
DIS paradigm, and are a direct generalization of the lat- 
ter. To access GPD's, it was proposed to study hard 
exclusive processes like deeply virtual Compton scatter- 
ing, (DVCS) e + p — >• e+p + 7 015] or meson electro- 
production, e + p — > e + p + p, bj, at high virtuality Q 2 
of the photon originating from the scattered lepton, and 
low momentum transfer t between recoiled and target 
nucleon. At present DVCS has been studied experimen- 
tally at HERA and Jefferson Lab [E2 QJ]. In- 
terpretability of hard exclusive processes in terms of the 
GPD's that are universal objects for all such reaction, is 
empowered by the collinear factorization theorem |15U16j 
that, similarly as for DIS, allows for a separation of the 
soft hadronic amplitude from perturbative, QCD process 
with the former leading to four GPD's. To the lowest 
order in the QCD coupling, a s , the full amplitude then 
corresponds to the handbag diagram depicted in Fig. [T] 
Paratactical applications, however, rest upon, the a priori 
unknown rate of convergence of the perturbative expan- 
sion. At low Bjorken-xs QCD corrections to the handbag 
diagram involve large logarithms in both a s log Q 2 and 
a s log 1/xb- While significant progress has been made 
in devising various resumation schemes [T7H25] , to date 
no first principle solution for the scattering amplitude 
exists. It is also accepted that the natural physical inter- 
pretation of the low- a: b DIS is quite different from that 
of parton model description of the valence region [21H5U] . 
That many orders in the a s expansion may been needed 
to describe the \ow-xb region is consistent with the ample 
evidence that in exclusive electroproduction nonpertur- 
bative phenomena play an important role in the nomi- 
nally perturbative domain. The structure functions at 
Iow-Xb have the behavior characteristic to Pomeron and 
Regge phenomena, while at fixed momentum transfer, 
exclusive photon or meson electroproduction cross sec- 




FIG. 1: Handbag diagram representation of the Compton am- 
plitude 

tions can be well fitted in terms of simple functions of 
Q 2 and the center of mass energy, W rather then Q 2 and 

x B [SH33I. 

Recently we have proposed a model in which the 
diffractive phenomena that are expected to govern the 
\ow-Xb DIS are incorporated at the parton nucleon 
level [3U |35]. As discussed above, at the QCD side, 
at low-xs resumation of gluon ladders leads to compli- 
cated evolutions equations. However since at large center 
of mass energy, hadronic amplitudes are known to have 
a universal Regge scaling, we employ this phenomena 
to construct an effective parton-nucleon amplitude. In 
terms of the QCD description of (THUS], in the model 
an infinite class of diagrams, i.e. those shown on the 
left panel in Fig. [2] is absorbed into the definition of the 
parton-nucleon blob and the resulting electroproduction 
amplitude is then computed from the handbag diagram. 
The model originates from a study of Regge phenomena 
at the parton level in the context of DIS [36l [37] . Such 
effective parton-nucleon amplitude gives the correct de- 
scription of low-x structure functions, surprisingly, how- 
ever, we have found that in the case of DVCS it breaks 
collinear factorization, i.e. Bjorken scaling while it nat- 
urally leads to the Regge-type scaling [Ml |35]. Upon 
closer examination, breaking of collinear approximation 
is not unexpected since it rests upon the assumption that 
parton-nucleon amplitude is a soft function of the invari- 
ant parton-nucleon energy, s. This is not the case if the 
amplitude has Regge-type, s", a > dependence on s. 
Such Regge-type scaling of exclusive amplitudes at large 
Q 2 and all ig as opposed to Bjorken-scaling was in fact 
predicted by Bjorken and Kogut in [55] , 

In this paper we focus on applicability of the the model 
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FIG. 2: Infinite class of perturbative gluon ladders (left) is 
expected to lead to Regge phenomena and are absorbed into 
the quark- nucleon amplitude (right). 



to DVCS in the HERA kinematics, Q 2 /W 2 « 1. For 
description of HERA data on DVCS at low-x two compet- 
ing formalisms are used, Regge models that operate with 
the soft and hard Pomeron trajectories, as for example 
in the color dipole and similar models piH^H |27| I3TH53] ■ 
and the GPD-based models. To be applied phenomeno- 
logically, the GPD-based models would include models 
for Regge-like background, see e.g. pMBX] , In general 
Regge background thus represents a systematic effect on 
the extraction of GPD's. Since both kinds of models are 
more or less successful in describing the HERA DVCS 
data, a question arises on whether the extraction of the 
GPD's is model independent. Moreover, if data allow for 
interpretation without GPD's, as in the model we study 
or the color dipole models, one may question the physical 
content of all these models. 

The paper is organized as follows. In the follow- 
ing section we discuss the DVCS amplitude in the 
handbag approximation and emerging properties of the 
parton-nucleon amplitude based on Regge phenomenol- 
ogy. Computation of the DIS an DVCS amplitudes is 
discussed in Sec. IIIII with more details included in the 
Appendix. Results and comparison with HERA data are 
presented in Se c. |IV| and followed by summary and con- 
clusions in Sec. [Vj 



II. COMPTON AMPLITUDE IN THE 
HANDBAG APPROXIMATION 

The hadronic Compton tensor is given by the matrix 
element of the time-ordered product of two electromag- 
netic currents, 



d 4 ze l 



'■{N\T\J v {z/2)J^{-z/2)]\N) 



(1) 



where q{q') is the four momentum of the incoming (out- 
going) photon. We will consider both the DIS pro- 
cess that corresponds to the forward virtual Compton 
scattering with both photons spacelike, q = q', q 2 — 
q 12 = -Q 2 < 0, and DVCS with q 2 < 0, q' 2 = 



and A = q — q' ^ 0. The currents are given by 
■^(z) = E,e 9 -W> W = Mz)l^ q {z) with ^ 
the quark field operator and e q the quark charge. Using 
the leading order operator product expansion we replace 
the product of the two currents by the product of two 
quark field operators and a free quark propagator be- 
tween the photon interaction points z/2 and —z/2, see 
Fig. [I] In this (handbag) approximation the hadronic 
Compton amplitude is then given by a convolution 

= ij 0^Q(K,q,6)A a0 (K,A,p,\,\') (2) 

of the quark Compton tensor 



Cp(K,q,S) 



(3) 
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a, ft being the Dirac indices, and the untruncated, with 
respect to the parton legs, parton-nucleon amplitude, 

= -iJ d 4 z e -^ z (p'A'|T[^ Q (z/2)^(-z/2)bA).(4) 

Following [331 [37] , we represent this amplitude as 

dfj, 2 



A a ,(K, A,p, A, A') - j {k , 2 _^ 2 + i£)(fc2 _^ 2+ . e) 
[(¥' + H)TK¥ + M)U S(PW u(p) (5) 



where Tf'? are constructed from Dirac 7-matrices and 



the available four-vectors p, A, k. The amplitude in 
gives the correct result in perturbation theory, 
e.g. for point-like quark-nucleon interaction. For partons 
bound inside the nucleon, however, A is expected to be 
suppressed at large- k 2 or k' 2 . This is achieved [M1I3T] . by 

so that in Eq.Q, 



1 



(fc' 2 - [i 1 + ie)(k 2 - [i 2 + it) 
1 



n {k' 2 -(j? +ie)(k 2 -fx 2 +ie)' 



(6) 



This method of softening the UV behavior guarantees 
current conservation. This would not be the case, for ex- 
ample, if the two propagators were absorbed into a soft 
quark-nucleon wave function. Furthermore, differentiat- 
ing the product of two propagators instead of differen- 
tiating each one separately ensures that the amplitude 
contains simple poles that enable to interpolate between 
the off- and on-shell quark-nucleon amplitudes. 
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FIG. 3: Direct and crossed contributions to the quark 
nudeon-scattering amplitude 



A. Quark-nucleon amplitude with Regge behavior 

We proceed by constructing the basis for the scatter- 
ing process N(p) + q(—k) —> N(p r ) + q(-k') shown in 



Figj3] We account for all possible Dirac-Lorentz struc- 
tures that can appear in four fermion operators. Fur- 
thermore we shall only consider those amplitudes which 
conserve the quark helicity since helicity-flip amplitudes 
are suppressed when integrated over in the handbag dia- 
gram by a power of p/W. The structures of interest thus 
involve ~ 7^,7^75 on the quark side only. Based on P, 
CP and CPT invariance, the quark-nucleon scattering 
amplitude can be decomposed in the basis of six inde- 
pendent tensors each then multiplied by a Lorentz scalar 
function, o», i = 1, . . . 6, 



A qN = qj a qN 



aila + a 2 - 



2M 



N ■ 



a3qi a i5qN'f a j 5 N + q-^q N 



a 4 + a 5 



M 



N ■ 



1 777T 9^7777^6, (7) 



2M 
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were the new four- vectors are defined by, K = (k + k')/2, 
P = (p + p')/2, and A = k' - k = p' - p. The amplitudes 
Oj are analytic functions of invariants s = (p — k) 2 — 
(P-K) 2 ,u = (p' + k') 2 = (P + K) 2 and t = A 2 , fixed 
by the condition s + u + t = 2M 2 + 2p? where p is the 
mass of the effective quarks (c.f. Eq.([5])) and we have 
explicitly put the quarks on the mass shell. The above 
basis is equivalent to the form used in [42] for elastic 
electron-proton scattering. In particular the amplitude 
multiplying 03 is chosen to be an axial vector but can be 
expressed in therms of qJPqNfiN used in pl2 . Moreover, 
fi. r y 5 /2M in front of a$ becomes proportional to 75 for 
on-shell particles, whereas f£/M multiplying 04 and 05 
reduces to \ijM. 

The scalar amplitudes a% have unitarity cuts in s and 
u and at fixed-i can be represented through a dispersion 
representation, 



tral functions, 



a,i(s,u,t,p 2 ) = (2tt) 4 



ds 



s — s 
subtractions 



(8) 



with the spectral function p^ being non-zero above some 
threshold values So(wo) m the respective channel. Next, 
we consider the phenomenological consequences of Regge 
exchanges on the asymptotics of the spectral functions 
at high s(u). For fixed t, we assume that the on-shell 
quark-nucleon helicity amplitudes follow Regge asymp- 
totics, i.e. they are proportional to s a ^ or u a W for 
large s or u respectively, a(t) being a Regge trajectory. 
Evaluating the asymptotic behavior of the amplitudes in 
Eq.Q and comparing with the expected behavior in the 
Regge limit we find the asymptotic behavior of the spec- 
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Note that in the pure collinear kinematics A M = 
(A + ,0,Oj_) (thus for A 2 = 0), and for massless quarks 
and proton, the matrix elements at 02,4,6 vanish identi- 
cally. Therefore, they generally have to be proportional 
to masses M, p or momentum transfer A 2 that is kept 
constant in Regge limit, and the above relations follow. 

An additional constraint on the behavior of the spec- 
tral functions comes from the Pomcranchuk theorem 
which implies that asymptotically s and u channel ampli- 
tudes become equal. The s — u crossing is implemented 
on the level of the quark-nucleon amplitudes according 
to 



K 
A 

7" 
7°7 5 



-> -K 

-> A 

-> C 7 Q C t = -7° 

-> C 7 Q 7 5 C t = + 7 a 7 5 , 



(10) 



with C denoting the charge transformation. For the spec- 
tral functions in Eq. |7| Pomeranchuk's theorem then im- 
plies, 



pf(s^ 00) 



-pf(«-> 00 ) 



for i 
for i 



3,4,6, 
1,2,5. (11) 
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We next introduce the C-even and C-odd combinations 
— ± p")/2 which asymptotically behave as, 



P 



Pi ~ s 

P2 

pi ~ s 
pt 



Ps 



s a2 

0-3 — 1 
,"5-1 



Pi 



Pi 

pi 

P3 
Pi 

Pi 

Pa 



(12) 



We notice that p~ and pf correspond to singlet (valence 
+ sea) and non-singlet (valence) GPD's, respectively. It 
is instructive to observe that according to Eq.(12), only 



singlet combinations may grow with s in the high energy 
regime, while the non-singlet ones necessarily vanish at 
high s. This fact, trivial in itself since it simply incor- 
porates the symmetry of the interaction of the nucleon 
with highly energetic quark and antiquark, has important 
consequence for collinear factorization. 

In Eq.Q, convergence of the dispersion integral at 
high energies is governed by asymptotic ener gy depen- 
dence of pi I s and p~ / 's 2 . Combining Eqs.(|12|),(|8|, it fol- 
lows that one can at most expect three subtraction con- 
stants, for 0,2, 04 and a§ [43 . The appearance of a finite 
subtraction constant that is energy- independent and thus 
has no exponential i-dependence would necessarily imply 
an appearance of fixed poles with very mild i-dependence 
in nucleon-nucleon and hadron-nucleon scattering. As it 
was noticed long ago [37J EI! ; the experimental data do 
not support such possibility and we will assume in the 
following that these subtraction constants are zero. 

The Pomeron can only contribute to the amplitude 
ax. The amplitude 03 has quantum numbers of an ax- 
ial vector ai-meson exchange which has the intercept 
a ai (0) « 0.5 and needs no subtraction. The amplitude 
05 is crossing-odd and needs no subtraction. 



and non-singlet GPD's. We combine Eqs.l]!]),^^ ^ 
to obtain the representation for the hadronic Compton 
amplitude 



i J dp 2 ds J d 4 Kj2u(p')T? 



u{p) 



s-(P- K) 2 



Tr 



s-(P 



Ky 

t 
2 



" \{K + A/2) 2 - p? + ie] [(K - A/2) 2 — p 2 + ie] 

Next, we will evaluate the contribution to the hadronic 
Compton amplitude from quark-nucleon amplitude pro- 
portional to ax, i.e. use u{p')Tf r u{p)Tj = u(p')j a u(p)j a 
(i = 1). This amplitude corresponds to Pomeron (and 
vector meson) exchange, so it should give the dom- 
inant contribution for DVCS at high energies where 
DVCS data from HI and ZEUS are available, We choose 
the kinematics [IS] as p^ = (p + ,0, 0±) and — 
(0, Q 2 /(2xbP + ),Q±), with the usual Bjorken variable 
xg = Q 2 /2pq. The trace in Eq.(13l can be evaluated 
using the collinear approximation 



w +p)t^(y+p)r 
1 



f(fcI+M 2 ) 



2x B P+* 



(K 



q+q \2 



ie {K- q Af) 2 



(14) 



III. REGGE EXCHANGE CONTRIBUTION TO 
DIS AND DVCS 

In this section, we will employ handbag formalism and 
relate the quark-nucleon spectral functions pi to singlet 



Note that the above trace calculation in collinear kine- 
matics is the same for forward (DIS) and non-forward 
(DVCS) case. Before we proceed, we notice that the trace 
in Eq.(14| is antisymmetric under exchanging K — > —K. 
This implies that only p~ spectral density contributes 
leading to, 



T£ = -Mg^^^u{p') 7 + u{p) I dp 2 ds I d±KI n 



x B 2P+ 



k\ + p 2 



x p x (s) 



1 



1 



[{K + A/2) 2 - p 2 + ie}[(K - A/2) 2 - p 2 + ie] 
1 1 



s - (P - K) 2 + ie s~(P + K) 2 + ie 



2+<L\2 



ie (K- 



(15) 



The fact that the above Compton amplitude depends on the collinear approximation: the positive C-parity of the 
the singlet spectral function p x only, is independent of Compton amplitude requires the C-even singlet combi- 
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nation p^ . On the contrary, the form factor, possessing We make the collinear approximation in the hard quark 
the odd C-parity only depends on the C-odd non-singlet propagators, 
combination p\ . 



A. DIS (7*p -> 7*p) 



We next evaluate the amplitude of Eq.(15) in the for- 
ward kinematics A = 0, q 2 = q' 2 = — Q 2 . 



T^A = 0) = -4< 

x f dp 2 ds J d A kI n 
1 



„.^ 2 1 

'" — u(p')j + u(p) 



x B 2p^ 
kj+p 2 
(fc 2 - p 2 + ie) 2 
1 



Pi ( s ) 



s — (p — k) 2 + ie s — (p + fc) 2 + ie 
1 1 



(k + q) 2 + ie (fc - q) 2 + ie 



1 



(fc + g) 2 +ie 
1 

(fc — g) 2 + ie 



x B /Q 2 



-Q 2 + ^ 



_g2 _ 

^ x B p+ 



fe+ 



^+ X B - le 



(17) 



(16) and obtain (we refer to the Appendix A for more details), 



T^(A = 0) = -4<gf — u{p')j + u{p)T{n) / dx(l - x) n+L dp\p z ) n dsp^ {s,t = 0, pf) 



2p+ 



dfc H 



K: - x B + ie K: + x B - ie 



[S(k+ -xp + )-8(k+ +xp+)] 



(n + 1 — x)p 2 + xs 
-(1 - x)p 2 - xs] n+1 

(18) 



To make a connection to the PDF's, we consider the 
imaginary part of this amplitude. Recalling that the 
imaginary part of forward Compton tensor proportional 
to — g^ v gives nWi — > 7r| J2 q e q[l( x ) ~ Q(~ x )]i we iden- 
tify the parton densities with integrals over the p~ or 
explicitly s and u spectral functions as, 



x B [q(x B ) - q(-x B )} = -8^ 2 r(n)(-l)" +1 (l - x B ) n+1 

[ d^d^ripK- 

I x B 
£+{n + l-x B )p 2 



x / dv 2 dap 2 ) n [p s 1 ( — ,0,p 2 )- P u 1 (- — Ap 2 



x B 



(£ + (l-ZB)/i 2 )" +1 ' 



(19) 



In the above, we changed the integration variable s to 
£ = x B s. Using the high energy asymptotics (c.f. Eq.Q) 
Pi' U ( s ) ~ s ap ~ 1 , with oip = 1 + e being the Pomeron 
trajectory, and pull the x B dependence out of the £- 
integral we obtain the experimentally observed asymp- 
totics F2(x B ) ~ x B ap ~ ^b 6 - This is the result for the 
singlet PDF. The non-singlet combination will depend on 
a similar integral with the non-singlet spectral function, 
which at high energy behaves as p + (s) ~ s ap ~ 2 , and cor- 



respondingly gives x B [q{x B ) + q(-x B )] 



x B , 



as expected. Evaluating the real part of the forward 
Compton amplitude we obtain the familiar result for DIS, 



1 



T^(A = 0) = — u(p') 7 +u(p) 



dx 



2x 

.x z — x B + ie 



(20) 



While the singlet PDF's at low x rise as x~ ap , the sin- 
gularity at x — > is cancelled by one power of x in the 
numerator of Eq.(20) which makes both the imaginary 



and real part of the integral finite [34] . 



B. DVCS (7*p — > 7f>): collinear approximation 



(p+ o,o_l), <f 



Next we evaluate Eq.(|15) in the DVCS kinematics, 
"(0,Q 2 /(2 ; r B p+),Q ± ) 



A M 



(— a:sp + , 0, 0j_), and choose now asymmetric integration 



variable fc, rather than K = ^ > 
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(21) 



Using the collinear approximation for the quark prop- 
agator exchanged between the two photons interaction 
points we obtain in the case of DVCS, 



1 
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x B /Q 2 



(k + qf 



-Q 2 



ie ^f-x B +ie 



(k-q 
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1 = x B /Q 2 

-k+ + ie -*i+ie" 
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(22) 



The DVCS amplitude in the collinear approximation is 
then given by, 



/ dx 
la 



1 



1 



x — xb + ie x — ie 



H + (x,xb), 



(23) 



and we refer the reader to the Appendix B for the details 
of the calculation. We identify the singlet GPD H(x,Xb) 
with, 

H+(x,x B ) = (I-xb/2) [ dy [ dz[q{z) - q{-z)} 

Jo Jo 

x 5(x — z — uxb(1 — z)) 
which satisfies the familiar normalization condition, 



(24) 



L 



1 f 1 

dxH + (x,XB) = (1 — xb/2) / dx[q(x) — q(—x)]. 

Jo 

(25) 



The factor (1— xb/2) in the definition of the GPD results 
from the prefactor 1/2P+ in the DVCS amplitude. Un- 
like DIS, in the presence of Regge asymptotics, the real 



part of the integral in Eq.(23) is divergent. This can be 



seen by first integrating the <5-function over x, and then 
performing the integral over y. In the limit z — > the 
real part of the integral 



dy 



1 



1 



z — xb + DXb(1 — z) + ie z + j/£g(l — z) — ie 

(26) 

is finite, and equal to ln(l — xb)/xb Then, given the 
Regge asymptotics of the PDF, [q(z) — q(—z)} ~ z~ ap 
the integral over z diverges. In the case of the DIS 
amplitude the quark propagator exchanged between the 
two photons in the sum of direct and crossed handbag 
diagram (c.f. Fi g. [T] ) leads to the factor of x in the 
numerator of Eq.(|20[). This does not happen in DVCS 
when one photon is soft and the sum of the two collinear 
propagators in the DVCS amplitude of Eq.(23) does not 
vanish when x — > and cannot compensate for the rise 
of the GPD at low x. We also note that in the case of the 
non-singlet GPD, the integral over x instead reduces to 
~ dxx 1 ~ ap and is therefore convergent. Thus conclude 
that for valence GPD's where Regge contributions are 
suppressed the collinear approximation is adequate 
and that part of the full DVCS amplitude would obey 
Bjorken scaling. As we show in the following section, 
inclusion of Regge contributions into singlet GPD's leads 
to Regge scaling. 



C. DVCS beyond the collinear approximation 

We will use the collinear approximation in the numera- 
tor only. We combine all four propagators together using 
Feynman parameters to obtain 



O 2 1 



= -8ig^ 



x B 2p+ 



x / d 4 k 



uy u J d/i I n I dsp^ (s)r(4) \ dxdydzil — x)(l — z) 

kl+p 2 



(27) 



([k + zq-(l- z)xp + y(l - x)(l - z)A} 2 - z(l - z)Q 2 (l - x/x B - y(l - x)) - (1 - z)[xs + (1 - x)n 2 ]) 4 

k\ +y 2 

([k - zq' - (1 - z)xp + y{\ - x)(l - z)A] 2 - z(l - z)Q 2 (x/x B + y(l - x)) - (1 - z)[xs + (1 - x)^ 2 }) 4 
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We report all the details of the algebra in the Appendix 
C, and quote here the final result, 



x 



1 



g + 3/i : x 
2 M 2 

"K + (i 



w(l - Zj?)|£ + M : 



- a; 



\t- 



g 2 + ^(i- 



y 2 + Q 2 - 



(28) 



* 

K + M 2 -^(1-^)]K + M 2 
(2- x B )(£ + ^ 2 +c^ B /2) 



2-x 



-u] 



B 



+ /i 2 + Q 2 - w] 



p 2 +Q 2 

(2-x B )(e + ^ 2 + Q 2 -^B/2) 

C + m 2 L[C + m 2 - w(i - ^s)] 2 [C + m 2 + H 2 K + m 2 + Q 2 + w (i - «b)] 2 K + m 2 + Q 2 - <A 2 



where we changed variables from s to £ = xs, from a; 
to w = Q 2 x/xb, and factored out the Regge asymp- 
totics of the spectral function as pj~(s) = s Q_1 /3 1 "(s) 
with j3 — ?> const, for s — > oo. Analyzing the above 
formula, we notice that integrals now converge. Impor- 
tantly, large values of ui do not contribute to the integral 
because of the explicit suppression factor (1 — xbuj/Q 2 ) 2 
and because of powers of lj in the denominator inside 
the bracket. The price to pay for this convergence is 
the appearance of the explicit scale dependence ~ p 2 in 
the expressions, as compared to the scale-independent re- 
sults obtained within the collinear approximation. This 
scale dependence is of no surprise since Regge behavior 
does introduce a scale. In the limit Q 2 /p 2 >> 1 it can 
be shown that the leading contribution of the Pomeron, 
a = ap to this integral is proportional to 



Tj 



1 



DVCS 
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X B 
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Q 2 



In the following Section, we will 
parametrization with the DVCS data. 



confront this 



IV. RESULTS AND COMPARISON WITH 
HERA DATA 



The result of the previous section, for Tjjvcs was ob- 
tained in the limit Q 2 — > oo at finite Q 2 the amplitude is 
finite but would require knowledge of the spectral decom- 
position of the quark-nucleon amplitude at finite ener- 
gies. When comparing to the experimental data at finite 
Q 2 we thus replace 1/Q 2 by a ~ 1/(1+Q 2 /<5q) with some 
characteristic scale Qq that we will determine from a fit. 
This is in accord with the experimental observation^ 



o~dvcs 



2\2ct-2 



2\t> 



(Q ) 



(30) 



with S rj —1.5 rather than —2. It is also this form that is 
used to describe data within phenomenological Regge (or 
color dipole picture-motivated) models [31] [33] . We will 



fit the HERA data using the following parametrization 
for the cross section 



00 



w 



i + Q 2 /Qi 



(31) 



(29) in Figs 



with Wo = 20 GeV. It is worth noting that using the reg- 
gized parton-nucleon amplitude in the handbag model we 
have effectively "derived" the parametrization proposed 
in [31]). 

We perform two fits. One is a combined fit to both HI 
[TO] and ZEUS QTJ 112] data. It gives a = 28 ± 4 nb, 
Qo = 1.51 ± 0.05GeV and a - 1 = 0.43 ± 0.03 and is 
shown in Figs IV with x 2 /d.o.f. = 2.01. The other, 
is an independent fit to HI and ZEUS data. For the fit 
to the HI data alone we obtain ao = 17 ± 3 nb, Qo = 
1.83 ± O.lGeV and a - 1 = 0.34 ± 0.05 and it is shown 
IV, with x 2 /d.o.f . = 1.2. For an independent 
fit to the ZEUS data alone we find ao = 41 ± 7 nb, Qo = 
1.49±0.06GeV and a-1 = 0.34±0.03 and it is shown in 
Figs|WJ[rVj with x 2 /d.o.f. = 1.1. We observe that both 



data sets are fitted well with the Regge form of Eq.(31 1, 



as it was found previously in color dipole or Regge based 
studies [3T]. However, the two data sets exhibit different 
normalization (the values of ao). As a result, performing 
a combined analysis we obtain a higher intercept. 



V. SUMMARY 

We presented an analysis of quark-nucleon scattering 
amplitudes. We considered a basis of six independent 
Dirac-Lorentz structures and discussed their Regge be- 
havior. In particular we have shown that the C-odd 
combinations of the direct and crossed channels (referred 
to as non-singlet combinations) follow different Regge 
asymptotics, as compared to the C-even (singlet) ones. 
Once embedded into the handbag diagram to describe 
the DVCS amplitude in hard kinematics, we show that 
only singlet combinations contribute, whereas the valence 



ZEUS W=I04 
ZEL'S W=S9 
H! W=82 
HI W=75 









FIG. 4: DVCS cross section as a function of photon virtuality, 
Q 2 for various cm. energies W (in GeV). In the upper panel, 
we confront the combined fit to the HI and ZEUS data. Solid 
lines a result of a fit to the combined ZEUS and HI data 
including both Q 2 and W dependence. The middle panel 
displays a similar fit to HI data alone, whereas the fits to 
ZEUS data alone are shown in the lower panel. 



combinations do not appear and require no a priori un- 
known subtractions. 

We focused on the contribution of a single Pomeron 
trajectory that dominates at high energies, and have 
demonstrate that while for DIS the handbag formal- 
ism leads to the known result, F%(xb) ~ x b" P : m the 
case of DVCS, the mismatch between quark propaga- 
tors leads to divergent integrals in the collinear approx- 
imation. If collinear approximation is not used, the 
model naturally leads to Regge-scaling for DVCS [38] 
with T DVCS ~ Q (2ap ~ 2) /x% p , with a P = 1 + e be- 
ing the Pomeron trajectory. Thus we have reproduce 




W [GeV] 

FIG. 5: W-dependence of the DVCS cross section for different 
values of Q 2 . The upper panel displays the comparison of the 
HI data to the combined fit to both data sets, whereas the 
second panel from top shows the ZEUS data vs. the same fit. 
The two lower panels confront individual fits to HI (second 
lowest panel) and ZEUS (lowest panel) to the corresponding 
data sets. 
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the form that phenomenological Regge models use to de- 
scribe DVCS, and we have illustrated its applicability by 
fitting the data from HERA. In he future we plan to ex- 
tend our phenomenological analysis to larger values of 
Bjorken xb, where DVCS was measured at Jefferson Lab 
[13], [14] . Since the JLab data is taken at much lower 
energies, however, the Pomeron trajectory alone is not 
expected to be sufficient and other trajectories will have 
to be studied. 
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Here we evaluate the forward Compton amplitude of 
Eq.pl, 



Tg(A = 0) = -Aig 
x p{ 



^ 2 1 

1 X B 2p* 
1 



(p'h + u(p) / dp 2 ds / d 4 kl 



1 



s — (p — k) 2 + ie s — (p + k) 2 + ie 



fcj + p 2 
(k 2 - p 2 + ie) 2 
1 



1 



(k + q) 2 +ie (k - q) 2 + ie 



(Al) 



Using the collinear quark propagators from Eq. (17) and introducing the Feynman parameter x, we obtain, 

1 1 



Tg = UgT^+uip^+uip) J dp 2 dsp^(s,A 2 ,p 2 ) J dk + dkrd 2 k ± 



K - x B +ie + x B - ie 



x {kl + ^W^TW 



i 



(p — k) 2 — s — ie (p + k) 2 — s — ie 



y_L 2 p^ 



x l dk 



w(p')7 + u(p)r(n + 3) / dx(l-x) n+L / dp 2 (p 2 ) n dsp^(s,A 2 ,p 2 ) 



^ - x B + ie + x B - ie 



x dk d 2 k± 



{kl+p 2 ) 



(kl+p 2 ) 



[(k - yp) 2 - ys - (1 - y)p 2 ] n + 3 [(k + yp) 2 - ys - (1 - y)p 2 } n+3 



(A2) 



Finally, Eq.(18) is obtained from Eq. (A2) after integrat 



ing over k , k± using, 
dk~d 2 k± 



dk d k± 



(k 2 +a 2 ) a 
(k 2 + a 2 ) a 



-ITT 



,T(a-2) 5(k+) 
T(a) (a 2 )"- 2 
2 r(a-3) 6{k+) 



T(a) (a 2 ) 



2W-3 ' 



(A3) 



The expression in Eq.(18) follows from Eq.(A2) after in- 
tegrating over k + . 

Appendix B: DVCS in collinear approximation 



We evaluate Eq.(15) in the DVCS kinematics, 
p» = (p+,0,0x), q» = (0,Q 2 /(2xbp+),Qi_)^ A^ = 
(— xbP + , 0, 0j_), and use k as the integration variable in- 
stead of K = (fr+fcO/2, 



ai 



^Q 2 i 



x B 2p~* 
1 



u{p)-/ + u(p) / dp? da / cTfc(fci + pf)I, 



1 



Pi 



s — (p — k) 2 + ie s — (p + k + A) 2 + ie 



[k 2 - p 2 + ie][(k + A) 2 - p 2 + ie] 
1 1 



(k + q) 2 +ie {k - q') 2 + ie 



(Bl) 
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We use the collinear approximation of Eq. (22) and com- to obtain 



bine the two quark propagators from the untruncated, 
quark-nucleon amplitude introducing an integral over a 
Fcynman parameter, 



1 



[k 2 - pi 2 + ie] [(k + A) 2 - fi 2 + ie] 



dy 



[(k + yA) 2 - fi 2 + ie} 2 ' 



(B2) 



K: = Wr^(p'h + u(p) J q dy I <I,M V ,] O-.A-./r) / <lk + <lk <fk_ 



2p 

x (ki + n 2 )i n 
1 



^-x B +ie ~^+ie 



1 



[{k + yA) 2 - /i 2 + ie] 2 
l 



(p — k) 2 — s — ie (p + k + A) 2 — s — ie 



Aig^—u(p'h + u(p) J dyT(n + 3) J dx(l - / dfi (fi ) n dsp^[ (s, A , /i ) 



x B + ie 



k+ 



x / dk d 2 k± 



(kl+^ 2 ) 



(B3) 



[(k -xp + y(l - x)A) 2 -xs-(l- x)p 2 ] n + 3 [(k + xp' + y(l - x)A) 2 - xs - (1 - x)p 2 ] n + 3 
Integrating over k~ , k± results in 

1 _, ,, _l , , f 1 , „ v„_i_i /" , o / . _, .o o, (n + 1 - x)p 2 + xs 



= -At:' 



2 9T^{p'h + <P) J q dyT(n) J dx(l-x) n+1 J d M 2 ( M 2 )"d S pr(s,A 2 V)p^ 



x)pL 2 — xs + ie] n+1 



dk+ 



1 



k+ 



x B + ie 



k+ 



[5(k+ -(x + yx B (l - x))p + ) - 8{k+ + (x(l - x B ) ~ yx B (l - x))p + )] 

(B4) 



The argument of the second (5-function can be brought tegration variables y — > 1 — y and k + 
to the same form of the first (5-function by changing in- Finally, the result reads 



T% = -4tt 



x 2 



(B5) 



(n + 1 — x)fi 2 + xs 



x B + yx B (l - x) + ie x + yx B {l-y) - it 



which corresponds to Eq. ( 23 1 with H defined in Eq. ( 24 ) . parameters 



Appendix C: DVCS beyond the collinear 
approximation 



We use the collinear approximation in numerator of 
Eq.( 13 1 and combine all four propagators using Feynman 
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x / d 4 k 



9l± 

x B 2p+ 



I 



wy + u I dp 2 I n J dsp 1 (s)T(4) j dxdydz(l — x)(l — z) 

k\ + a 2 



(CI) 



([k + zq-(l- z)xp + y(l - x)(l - z)A} 2 - z(l - z)Q 2 (l - x/x B - y(l - x)) - (1 - z)[xs + (1 - x) M 2 ]) 4 

kj+p 2 

' ([k - zq> - (1 - z)xp + y(l - x)(l - z)A} 2 - z(l - z)Q 2 (x/x B + y(l - x)) - (1 - z)[xs + (1 - x)/z 2 ]) 4 



Integration over d*k results 



o 2 uv Q 1 

8tt g 1 ^ — ^tx u T 



+ »/' 



d/j 2 /i 4 I n -2 I dsp 



Jo 



dxdydz{\ — xf 



x B 2p+ 

l-z l-z 
[xs + (1 - x) M 2 + zQ 2 (l - x/x B - yjl - x))] 3 ~ [xs + (1 - x)y? + zQ 2 (x/x B + y(l - x))] 3 

3( M 2 + z 2 Q 2 ) 3(^ 2 +z 2 Q 2 ) 

[xs + (1 - x)fi 2 + zQ 2 (l - x/x B - y(l - x))] 3 [xs + (1 - x)^ 2 + zQ 2 {x/x B + y(l - x))} 3 



(C2) 



Next the y integral can be done to obtain, 
1 1 



8n 2 q^ v — - — - wy 
y± x B 2p+ 



+u J dp 2 u A I n _ 2 J dspi(s)T(3) dx{l-xf 



1 



dz 



1 



+ 



[xs + (1 - x)fi 2 + zQ 2 (x - x/x B )] 2 [xs + (1 - x)p 2 + zQ 2 (l - x/x B )) 2 
1 1 



[xs + (1 - x)p? + zQ 2 (l -x + x/x B )] 2 [xs + (1 - x) f i 2 + zQ 2 x/x B } 2 



+2(p 2 +z 2 Q 2 ) 
+ ■ 



1 



1 



[xs + (1 - x)p 2 + zQ 2 (x - x/x B )] 3 [xs + (1 - x)p 2 + zQ 2 (l - x/x B )} 3 
1 1 



[xs + (1 - x)n 2 + zQ 2 (l -x + x/x B )} 3 [xs + (1 - x)p? + zQ 2 x/x B } 3 
and finally z integral yields, 



+ 



£+(3-x)/x 2 

v 

[£+(l-aO/x 



In 



212 



g + (1 - x)p 2 + Q 2 (x - x/ijj) £ + (1 - x)p 2 + zQ 2 (l - x + x/x B ) 
£ + (1 - x)p 2 + Q 2 x/x B £ + (l-x)n 2 + zQ 2 (l-x/x B ) 
1 1 



+ 



+ 



£ + (1 - x)p 2 + Q 2 (x - x/xb) £ + (1 - x)p 2 + Q 2 (l - x/x B ) 

1 



£ + (1 - x) M 2 + Q 2 (l - a: + x/x B ) £ + (1 - x)^ 2 + Q 2 x/x B 
/i 2 + Q 2 



£ + (i-x) M 2 



1 



1 



+ 



[£ + (1 - x) M 2 + Q 2 (x - x/a; B )] 2 [£ + (1 - x) M 2 + Q 2 (l - x/xb)] 2 
1 1 



[£ + (1 - x)^ 2 + Q 2 (l - x + x/xb)] 2 K + (1 - ^)M 2 + Q 2 x/x B } 2 

I 



(C3) 



(C4) 



where we changed variables from s to £ = xs and fac- since <~ £ a_1 and the expression in the curly bracket 

tored out the Regge asymptotics of the spectral function drops at least as l/£ 3 . Instead, the x integral is peaked 

as Pi (s) = s" -1 /3i (s) with j3 — > const, for s — > oo. To at x — > 0, and we can therefore neglect x in terms propor- 

proceed, we observe that the integral over £ is convergent tional to (1 — x). The divergent behavior of this integral 
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obtained in collinear approximation for the propagators gral of the type J Q dxx~ a . To ensure convergence, we do 

can obtained the formal limit Q 2 — > oo Then, the expres- not make this approximation. Changing finally t he in te- 

sion in the curly bracket becomes Q 2 -indcpcndcnt, and gration variable itow = Q 2 x/xb, wc obtain Eq.(28). 
proportional to ~ ln(l — xb) leading to a divergent inte- 
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